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SUMMARY 

This  research  program  focusses  on  the  electronic  structure  aspects  of  radiationless  decay 
processes  related  to  the  stability  and  formation  of  high  energy  density  materials.  We  have  also 
begun  theoretical  studies  on  electronic  nonadiabaticity  in  ion-molecule  reactions  These  studies  are 
enabled  by  a  unique  system  of  electronic  structure  algorithms,  the  BROOKLYN  programs,  which 
we  have  developed  over  the  last  decade.  These  programs  piuvide  advanced  capabilities  for  the 
study  of  the  electronic  structure  aspects  of  spin-forbidden  and  spin-allowed  electronically 
nonadiabatic  processes.  The  methodology  we  have  developed  is  based  exclusively  on  large  scale 
configuration  state  function  expansions  (lO^-lO^  terms).  These  methods,  which  we  believe 
define  the  state  of  the  art  in  these  areas,  have  permitted  us  to  make  significant  contributions  to  the 
understanding  of  electronically  nonadiabatic  processes.  Problems  of  particular  relevance  to  the 
high  energy  density  materials  program  include  studies  of  the  stability  of  the  potential  energetic 
species  a-N202,  tetrahedral  N4  and  the  dication  HS2+.  We  have  also  initiated  a  research  program 
in  ion-molecule  chemistry  reporting  for  the  first  time  a  seam  of  crossings  of  two  states  of  the  same 
symmetry  for  the  prototypical  charge  transfer  reaction  +  NO  ->  H  +  NO+. 

Finally  a  new  phase  of  program  development  has  been  initiated  with  the  goal  of  extending 
our  capabilities  for  the  direct  determination  of  actual/allowed/avoided  crossing  hyperlines  and 
hyperplanes.  The  algorithms  under  development  will  permit  systematic  determination  of  (i) 
actual/allowed  crossing  seams  passing  through  a  minimum  energy  crossing  point  and  (ii)  avoided 
crossing  seams  as  a  function  of  an  arbitrary  set  of  internal  coordinates.  This  program 
development,  which  will  be  completed  during  our  present  AFOSR  sponsored  research  program, 
will  provide  powerful  new  tools  for  the  study  of  spin-forbidden  and  spin-allowed  electronically 
nonadiabatic  processes  in  polyatomic  systems. 
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TECHNICAL  REPORT 

One  of  the  principal  aspects  of  our  research  program  has  been  the  development  of  new 
theoretical  tools  for  treating  spin-forbidden  and  spin-allowed  electronically  nonadiabatic 
processes.  These  tools  have  then  been  used  to  study,  the  stability  and/or  decay  mechanisms  of 
potential  high  energy  density  materials(HEDM’s)  and,  ion-molecule  chemistry.  We  believe  that  as 
a  result  of  computational  limitations  the  electronic  structure  aspects  of  these  important  classes  of 
chemical  problems  had  received  inadequate  treatment.  The  methodology  we  developed  is  based 
exclusively  on  large  scale  configuration  state  function  expansions  (10^-10^  tenns).  These 
methods,  which  we  believe  define  the  state  of  the  art  in  these  areas,  have  permitted  us  to  make 
significant  contributions  to  the  understanding  of  radiationless  decay  processes  and  electronically 
nonadiabatic  chemical  reactions.  Below  we  outline  the  methods  we  have  developed  and  their 
applications  in  the  HEDM  and  ion-molecule  chemistry  programs. 

I:  METHODS:  SPIN-FORBIDDEN  PROCESSES 

Our  studies  in  the  area  of  spin-forbidden  chemistry  employ  the  Breit-Pauli  approximation  ’ 
and  consider  those  interactions  which  lift  the  2S+1  degeneracy  associated  with  the  total  electron 
spin  quantum  number,  S.  These  interactions  are  the  full  microscopic  spin-orbit  interaction  H^; 

1=1  ■./=! 

and  the  dipolar  spin-spin  interaction  H^s 

=  1.2 
■./=l 

where  h^(i),  hsoo(i,j),  and  h*s(i,j)  are  the  microscopic  spin-orbit,  .spin-other-orbit,  and  dipolar 
spin-spin  interactions*  and  we  write 
(A)  Methodology 

Our  treatment  of  these  interactions  is  unique  in  two  ways.  We  use  the  symbolic  matrix 
element  methods  originally  introduced  in  the  treatment  of  the  nonrelativistic  Bom  Oppenheimer 
hamiltonian  H®  by  Liu  and  Yoshimine^  to  evaluate  the  matrix  elements  of  H^^and  in  the  N- 
electron  configu.ation  state  function  (CSF)  space.  This  enables  the  use  of  large  scale  CSF  spaces 


3 


(10^-10^  terms)  to  characterize  these  interactions.  These  spaces  are  one  to  two  orders  of 

magnitude  larger  than  those  accessible  using  previously  existing  methods  ^  and  this  capability 

has  improved  significantly  the  reliability  of  theoretical  work  in  this  area. 

The  second  unique  aspect  of  our  approach  involves  the  manner  in  which  the  relativistic 

wavefunction  is  obtained  from  the  zeroth  order  nonrelativistic  (configuration  interaction  (Cl)) 

7 

wavefunction  '  Within  the  context  of  the  Breit-Pauli  approximation  the  relativistic 
wavefunction,'F(I),  is  obtained  from  the  zeroth  order  nonrelativistic  wavefunction,  4^0(1),  using 
perturbation  theory.  In  order  to  allow  for  the  possibility  of  intersections  or  near  intersections  of 
potential  energy  surfaces  we  have  developed  an  approach  based  on  quasi-degenerate  perturbation 
theory.  The  description  which  emerges  from  this  formalism  expresses  the  electronic 
wavefunctions  as  a  linear  combination  of  dressed  diabatic  wavefunctions  f'P‘^(I)] 

=  +  1.3 

where  'PO(K)  satisfies  the  nonrelativistic  Schrodinger  equation 

[HO-EO(I)]'PO(i)  =  o  1.4 

and  the  dressing  term,  '¥^^(1)  represents  the  perturbation  of  the  (diabatic)  quasidegenerate  state 
'{'0(1)  by  states  outside  the  quasidegenerate  space  and  satisfies 

[//”  -  £'’(/)]'{'''“(/)  =  1 .5 

where  Q  projects  onto  the  orthogonal  complement  of  the  quasidegenerate  space.  The 
quasidegenerate  dressed  diabatic  states  '{'<1(1),  the  analogues  of  the  'POfl)  in  eq.  1 .4,  are  mixed  by 
the  Breit-Pauli  interaction  H®P.  This  mixing  is  given  by  a  secular  equation  with  the  dimension  of 
the  quasidegenerate  space  and  having  the  form: 


'Hu-E 

... 

"O' 

0 

• 

"2.2-^ 

= 

.0. 

where  H,  j  =  +  //®'’j'P‘'(y)J  .  Since  the  coupling  between  the  states  'P‘1(I)  originates 

from  a  potential  term  rather  than  the  derivative  coupling  terms,  to  be  discussed  sub.sequently, 
the  'P<1(I)  are  referred  to  as  diabatic  states. 
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This  approach  is  to  distinquished  from  the  commonly  used  approach  ^  in  which  eq.  1.5  is 
replaced  with  a  spectral  expansion  of  'Pq‘'(/) 


£'‘{K)-E\l) 


1.7 


This  procedure  is  unsatisfactory  owing  to  the  need  to  drastically  truncate  the  sum  in  eq.  1 .7.  The 
use  of  eq.  1.5  in  lieu  of  eq.  1.7  has  improved  considerably  our  ability  to  predict  rates  for  spin- 
forbidden  dipole-allowed  radiative  decay  processes. 


n:  NONADIABATIC  PROCESSES 

A:  Determination  of  Seams  of  Avoided/ Actual  Crossings  Seams  Using  Analytic  Derivative 

Techniques 

When  considering  the  stability  of  energetic  local  extrema  on  the  ground  potential  energy 
siuface  one  must  consider  the  possibility  of  interactions  with  nominally  excited  potential  energy 
surfaces  as  these  interactions  may  lead  to  dissociation  channels  for  the  metastable  energetic 
species.  Yet  in  the  absence  of  simplifications  attributable  spacial  symmetry  such  regions  of 
nuclear  coordinate  space  are  difficult  to  isolate.  On  the  basis  of  these  simple  observations  rests  the 
need  for  an  algorithm  for  the  determination  of  regions  of  close  approach  of  two  potential  energy 
surfaces.  Regions  of  close  approach  can  be  classified  as  actual/allowed/avoided  surface 
crossings.  Actual  (allowed)  surface  crossings  represent  the  intersection  of  two  potential  energy 
surfaces  of  the  same  (different )  symmetry.  The  set  of  all  such  points  forms  a  surface  of 
dimension  N-2  (N-1)  where  N  is  the  number  of  internal  nuclear  degrees  of  freedom.  Actual 
crossing  are  permitted  by  the  multidimensional  breakdown  of  the  vonNeumann-Wigner 
noncrossing  rule.^^  Since  the  multidimensional  breakdown  of  the  noncrossing  rule  permits  but 
does  not  guarantee  the  existence  of  actual  crossings,  actual  crossing  seams  can  merge  smoothly 
into  avoided  crossing  seams  of  the  same  dimension.  Our  capability  for  the  determination  of  seams 
of  avoided/actual  surface  crossings  rests  on  an  isomorphism  connecting  the  equations  used  in 
constructing  energy  gradients,  energy  difference  gradients  and  the  (Cl  contribution)  to  the  first 
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derivative  nonadiabatic  coupling  matrix  elements.*  *  The  essential  elements  of  this  isomorphism 
are  outlined  below. 

We  arc  required  to  construct  a  generalized  first  derivative  with  respect  to  the  internal 
coordinate  Ra.  This  quantity  has  the  form 


q“(R)  =  QS(R)  +  Q|J°(R)  2-1 

where  the  first  term  gives  the  contribution  attributable  to  changes  in  the  atomic  orbital  integrals 

e;=lpj'A,“  +  2-2 

i.j  ij.i.l 

and  the  second  term  gives  the  contribution  attributable  to  changes  in  the  molecular  orbitals 

Gr=ZL“U“  .  2.3 

i.j 

Here  and  defined  further  below,  are  generalized  one  and  two  particle  density  matrices  for 
multireference  Cl  wavefunctions  and  'PO(j)  and  the  generalized  Lagrangian  is  given  by 

Lt  =2[Xp“  +  2£PS,(mjllkl)l  .  2.4 

.  j  i.*'.' 

In  the  above  h®  and  (ijllkl)oi  are  the  derivatives  of  the  standard  one  and  two  electron  integrals  with 
respect  to  Ra,  and  the  U°’(R)  represent  the  derivati  ves  of  the  molecular  orbital  coefficients,  that  is 


^^,(r;R)  =  5;^,('-iR)G,°+I  ^2:/('';R)  <,<R) 

'a  J  i 


where  i|>(r;R)  [X(r;R)]  are  the  molecular  (  atomic  )  orbitals  and  t(R)  expresses  the  molecular 
orbitals  in  terms  of  the  atomic  orbital  basis.  The  U®(R)  are  determined  from  the  coupled  perturbed 
state  averaged  MCSCF  equations.*^’ 

Eqs.  2. 1-2.4  constitute  the  aforementioned  isomorphism.  The  desired  derivative  quantities 
are  obtained  by  replacing  the  general  density  matrices  p  and  P  with  specific  density  matrices  as 
given  below: 

(a)  Energy  Gradient  ^-£'(R):  standard  density  matrices*^’ 

dR„ 

pH  tI'  and  Pli,  rj;^„ 
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(b)  Energy  Difference  Gradient  — — A£'^(R) :  difference  density  matrices 

oR„ 


andl^y-»ri|„-r” 


'U 

d 


ijkl  ^  ijkl 


where  AEU(R)=Ei(R)-Ej(R);  and 

(c)  Derivative  Couplings  /4'°(r;R)l-^*Py  (r;R)\  :  transition  density  matrices 


SR. 


p|'-»  yI’  and  ^  r”„. 


ijkl  ^  Mjkl’ 


q^\R)/AEU(R)-4 

with  the  complete  first  derivative  nonadiabatic  coupling  matrix  element  ^(R)= 
^'F;(r;R)|  J-'P]’(r;R)^  given  by;*^- 

C(R)  =  X'(^)  +"7i^R) 


where 


and 


CSF 


d“(R)  =  (^,(r;R)|^^/r;R) 


2.6 

2.7 

2.8 


The  antisymmetric  matrix  d“(R)  is  readily  determined  from  eq.  2.5  once  the  U“(R)  are  available. 

The  direct  determination  of  the  energy  difference  gradient  is  the  key  to  our  algorithm  for  the 
efficient  determination  of  actual/avoided  surface  crossings.  * '  This  intersurface  feature  is  obtained 
by  minimizing  AE^(R)2  using  a  Newton-Raphson  procedure.  The  requirement  that  AE^CR)^  be  a 
minimum  yields: 


4-A£"(R)'=2A£"(R)ga(R)  =  0  2.9 

Eq.  2.9  shows  that  both  avoided  and  actual  surface  crossings  can  be  obtained  from  the  same 
algorithm  and  evinces  the  role  of  the  energy  difference  gradient  in  such  a  procedure. 
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B:  Determination  of  the  Minimum  Energy  Crossing  Point  on  the  Surface  of  Intersection  of 

Two  States  of  Different  Spin-Multiplicity  Using  Analytic  Derivative  Techniques 

When  considering  a  nonadiabatic  processes  involving  states  of  different  spin-multiplicity 
one  is  only  concerned  with  allowed  surface  crossings.  The  minimum  energy  point  on  the  surface 
of  intersections  represents  a  key  bottleneck  along  the  minimum  energy  path  in  a  spin-forbidden 
process  and  frequently  represents  the  barrier  for  this  process.  Thus  the  direct  ''haracterization  of 
this  point  is  an  essential  first  step  in  any  r.eatment  of  a  spin-forbidden  electronically  nonadiabatic 
process. 

The  straightforward  procedure  that  has  been  used  in  the  past  to  locate  the  minimum  energy 

to 

point,  the  indirect  determination  of  this  point,  initially  determines  the  crossing  surface  and  then 
characterizes  its  minimum.  This  procedure  is  computationally  costly  since  for  a  system  with  N 
internal  degrees  of  freedom  a  crossing  subsurface,  or  hyperplane,  of  dimension  N-1  must  be 
determined  and  analyzed.  Thus  the  determination  of  this  point  represents  a  significant 
computational  bottleneck.  However  this  bottleneck  can  be  avoided.  The  minimum  energy 
crossing  point  can  be  determined  directly,  that  is  without  prior  determination  of  the  surf  ice  of 
intersection  itself,  using  a  Lagrange-Newton  procedure*^  which  minimizes  the  quantity 

LJ-I  (K,X)  =  E°(R)  -t-  MEj(R)  -  Ef(R)]  2. 10a 

=  Ej®(R) >J^iEJl(R)  2.10b 


where  R  denotes  the  nuclear  coordinates  and  A,  is  a  Lagrange  multiplier.  Expanding  eq.  2. 10  to 
second  order  yields  the  Lagrange-Newton  equations*  *'  *^' 


"’W(R,A) 


g"(R) 
0 


^^6R  ^ 


A 


V(R)+V'(R)' 

.  4£"(R)  , 


where  8R  =  R'-R,  dX  =  X'  -  X, 


w"(Ra)  = 


d^L,jiR,X) 

dR„9Rp 


J-(s;(R)+Ag"(R)] 


2.11 


2. .12 


and  the  gradients  are  defined  as 
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dA£„{R) 


and 


gim 


dR., 


2.13b 


We  have  implemented  an  algorithm  for  the  solution  of  eq.  2. 1 1  in  which  the  gradients,  cq. 
2.13,  are  evaluated,  for  MCSCF/CI  wavefunctions,  using  analytic  derivative  techniques  (eqs  2.1- 
2.4).  *  is  deteruiiucd  using  a  forward  or  a  centered  difference  of  +  ?>.gp^  that  is; 

W"(R.A)  =  {[4(R  +  £r)  +  Ag"(R  +  £l“)]-[4(R)  +  Ag;^(R)]}/f  2.14a 


or 

=  {[jJCR  +  d“)  +  A*;(R  +  £l“)]  -  [s;(R  -  £l“ )  +  As;(R  -  £l“  >1}  /  2f  2.  Mb 

where  1“  is  a  unit  vector  along  the  direction  Rq. 

The  use  of  analytic  gradient  techniques  permits  the  solution  of  eq.  2. 1 1  in  ihc  full  nuclear 
coordinate  space.  An  analysis  of  the  nuclear  configuration  and  wavefunctions  at  the  minimum 
energy  crossing  point  can  provide  valuable  mechanistic  insights  into  a  spin-forbidden  reaction." 
Ill.  APPLICATIONS 
A.  High  Energy  Density  Materials 

(i);  Stability  of  Dicadons  with  Respect  to  Spin-Forbidden  Radiadve  Decay 

In  a  recent  combined  theoretical/experimental  study  at  the  University  of  Colorado  at 
Boulder  Miller  et  al?^  reported  the  existence  of  the  dication  HS^'*'.  The  vibrational  levels  of  the 
.state  observed,  the  X^n  electronic  state,  must  be  meiastable  with  a  lifetime  greater  than  10  jisec  to 
be  detected  in  their  experiment.  Using  large  scale  complete  aaive  space(CAS)"''‘"^/configuration 
interaction(CI)  wavefunctions  these  authors  determined  the  potential  energy  curve  for  the  X^O 
state  (see  figure  below).  They  showed  that  the  barrier  in  the  X^O  potential  energy  curve, 
attributable  to  an  avoided  crossing  between  states  correlating  with  H‘*'S'^  and  HS^"^.  can 
support  eight  quasibound  vibrational  levels.  Tunneling  on  the  X^FI  potential  energy  curve  limits 
the  lifetime  of  the  vibrational  levels  with  v>5  to  lOOpsec  or  less.  However  the.se  authors  also 
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determined  that  a  dissociative  electronic  state  crosses  the  bound  X^n  state  at  R=  1 .8 1 A 
(3.42ao)  at  an  energy  above  that  of  X^n,  v=l  level.  This  crossing  provides  an  additional  decay 
channel  for  the  vibrational  levels  of  the  X^fl  state  resulting  from  the  spin-orbit  perturbation 
X2n~l4l-. 

As  part  of  collaboration  with  the  Boulder  group  we  undenook  a  study  of  the  efficacy  of  this 
mechanism  using  a  combination  of  ab  initio  electronic  structure,  and  coupled  state  nuclear 
dynamics,  calculations.  The  following  figure  summarizes  the  results  of  our  ah  initio  calculations 
which  determined  both  the  potential  energy  curves  for  the  X^fl  and  HL-  states  as  well  as  the 
requisite  spin-orbit  coupling  l^InX-n^)  s  ('F(l^In)  I . 
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10000 
5000 
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-5000 
-10000 
-15000 
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Using  this  data  the  radiationless  decay  lifetimes  were  determined  within  a  Hund's  case  (a) 
approximation.  The  most  important  conclusion  of  this  study  is  that  the  1^L~  ~  X^fl  state 
interaction  significantly  shortens  the  single  state  tunneling  lifetimes  of  the  vibrational  levels  of 
the  X^n  state  reported  by  Miller  et  al.  Despite  the  fact  that  the  radiationless  decay  of  the  v=0.  1 
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levels  is  due  to  spin-orbit  induced  tunneling  we  find  that  only  the  v=0  level  will  he  obsen  able 
on  the  time  scale  of  the  above  noted  experiment. 

(ii) ;  Spin-forbidden  Radiationless  Decay  in  Tetrahedral  N4. 

The  potential  energetic  material  tetrahedral  N4  has  been  the  subject  of  several  recent  high 
quality  theoretical  studies^^’  that  have  demonstrated  that  this  species  is  in  fact  a  local 
iTunimum  on  the  lowest  ^A'  potential  energy  surface.  However  this  species  has  yet  to  be 
synthesized.  The  barrier  to  spin-allowed  decay  of  N4 

N4(1Ai)  N2(XlXp  +  N2(XlXp  3.1a 

is  approximately  61kcal/mol.^^  It  has  been  suggested^  that  the  spin-forbidden  decay  channel 

N4(‘ Ai)  N4(3A")  -4  N2(XlIp  +  N2(A5I;)  3. 1  b 

limits  the  liletime  of  tetrahedral  N4  and  is  therefore  responsible  for  the  inability  to  synthesize  this 
species. 

The  decay  channel  given  by  eq.  3.1b  had  not  been  previously  investigated.  Our  treatment 
of  this  system^^  found,  using  eq.  2.11  and  multireference  Cl  expansions  with  over  400,000  terms, 
that  the  minimum  energy  crossing  point  of  the  ^A'  and  ^A"  potential  energy  surfaces  is  only 
28kcal/mol  above  tlie  ^  Ai  minimum  on  the  ^  A'  potential  energy  surface  corresponding  to 
tetrahedral  N4!  This  important  finding  means  that,  (i)  the  lowest  vibrational  levels  of  N4  are 
unlikely  to  be  predissociated  by  this  mechanism,  and  (ii)  spin-forbidden  radiative  decay  must  be 
considered  in  any  realistic  treatment  of  the  lifetime  of  tetrahedral  N4.  We  intend  to  pursue  this 
point  further  using  a  golden  rule  wavepacket  method  for  determining  the  radiationless  decay  rate.'^* 
This  method  should  be  computationally  tractable  for  a  polyatomic  system  since  the  wavepacket 
n^  only  be  propagated  for  a  shon  time  on  the  dissociative  surface.^' 

(iii) :  On  the  Stability  of  a-N202 

The  asymmetric  dimer  of  NO,  a-N202,  had  previously  been  suggested  by  Michels  and 
Montgomery  as  a  potential  energetic  material.^^  In  collaboration  with  Mark  Gordon's  group  at 
Iowa  State  University  and  Michels  and  Montgomery  the  stability  of  this  system  with  respect  to  the 
exoergic  but  spin-forbidden  radiationless  decay  process 
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a-N202(*A’)  ->  N202(3A")  N20(Xll+)  +  0(3p)  3,2 

was  considered.  Large  scale  multireference  configuration  interaction  wavef unctions,  comprised 
of  approximately  300,000  -1,400,000  configuration  state  functions,  based  on  double  zeta 
polarization  and  triplet  zeta  polarization  bases  were  used  to  study  this  process.  The  minimum 
energy  crossing  of  the  ground  singlet  *A’,  and  the  lowest  excited  triplet  ^A",  potential  energy 
surfaces  was  determined  as  was  the  interstate  spin-orbit  coupling. 

a-N202  is  perhaps  best  described  as  0(*D)  datively  bound  to  N20(XiZ"*').  A  low-lying 
triplet  state  correlating  with  N20(XiX+)  -f  0(3P)  was  found  to  intersect  (at  a  geometrical 
arrangement  labelled  MEX)  the  ^  A'  potential  energy  surface  only  1 .0  -  2.0  kcal/mol  above  the  a- 
N2O2  equilibrium  structure,  denoted  MIN(*A’).  MIN(iA')  and  MEX  differ  only  in  R(0*-0“) 
which  increases  by  O.I2A  going  from  MIN(lA’)  to  MEX.  The  spin-orbit  coupling  //"'  was  found 
to  be  large  in  this  region,  H(°  ~  75cm-f  The  structural  similarity  among  MIN(*A'),  MEX  and 
isolated  N20(XlZ'‘‘)  permits  the  spin-forbidden  predissociation  to  be  studied  using  a  simple  one¬ 
dimensional  model.  This  model  predicts  that  a-N202  will  be  rapidly  (on  the  order  of  picoseconds) 

dissociated  to  N20(X^Z+)  and  0(3p)  as  in  eq.  3.2  and  thus  is  NOT  a  HEDM  candidate. 

B.  Ion-Molecule  Chemistry 

The  charge  transfer  reaction, B+  +  A  B  A+,  is  perhaps  the  archetypical  example 

-ic 

of  an  electronically  nonadiabatic  process  involving  a  transition  between  the  electronic  potential 
energy  surfaces  corresponding  to  B'*'  -t-  A  and  B  +  A'*’.  Such  processes  play  an  important  role  in 
the  chemistry  of  the  upper  atmosphere. 

(i):  Actual/Allowed  Crossing  Seams  in  the  Reaction  H'*'  -1-  NO  ->  H  +  NO"*" 

The  location  of  the  nuclear  configurations  for  which  the  potential  energy  surfaces  are  in 
close  proximity  is  of  prime  importance  in  elucidating  the  mechanism  of  a  charge  transfer  process. 

In  the  frequently  studied  charge  exchange  systems  (ArH2)'''  and  (ArN2)'^  (Ref  36 )  the  regions  of 
close  approach  are  characterized  by  symmetry  allowed  surface  crossings  of  the  relevant  potential 
energy  surfaces.  Here  symmetry  allowed,  or  just  allowed,  surface  crossings  refer  to  crossings  of 
adiabatic  electronic  states  of  different  symmetry.  Crossings  of  potential  energy  surfaces 
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corresponding  to  adiabatic  electronic  states  of  the  same  symmetry  are  normally  considered  rare  and 
are  usually  ignored  in  discussions  of  charge  transfer  phenomena.  Actual  crossing  seams  are 
difficult  to  locate  because  they  represent  surfaces  of  intersection  of  dimension  N-2  whereas  for 
allowed  crossings  the  surfaces  of  intersection  are  of  dimension  N-1.  In  fact  it  was  not  until  1990 
that  instances  of  actual  crossing  scams  based  on  MCSCF  (for  03)^^  and  MCSCF/CI  (for  He- 
H2)^^  wavefunctions  were  reported.  These  examples  had  been  the  only  reported  instances  of 
actual  crossing  seams  based  on  MCSCF/CI  wavefunctions.  No  examples  of  actual  crossing  .seams 
relevant  to  charge  transfer  reactions  existed  despite  the  long  history  of  theoretical  studies  of  this 
class  of  processes.^^ 

In  our  initial  study  of  the  portotypical  charge  transfer  system ,  HNO*’ ,  an  actual  crossing 
seam  of  two  potential  energy  surfaces  of  the  same  symmetry  relevant  to  the  charge  transfer  reaction 
H+  +  NO(X2n)  H(2S)  +  NO+(XlI+)  was  located.^^  This  crossing  seam.  R(R)s(R.r(R). 
7(R)]  occurs  for  general  Cg  geometries.  It  represents  the  intersection  of  the  1,2  2  A'  potential 
energy  surfaces  which  correlate  asymptotically  with  H(2S)  +  NO‘*'(X^I''')  and  H+  +  NOfX^ri) 
respectively.  For  R~2.45-3.65ao  ihe  seam  is  exoergic  relative  to  the  reactant  channel  asymptote, 
H+  +  NO(X2n).  The  present  results  complement  a  previous*^^  determination  of  two  allowed 
crossing  seams  in  Coov  symmetry.  Thus  these  results  provide,  for  the  first  time,  an  example  of  (i) 
an  actual  (  as  opposed  to  an  allowed)  crossing  seam  for  a  charge  transfer  reaction  and  (ii)  a  system 
for  which  both  allowed  and  actual  -  energetically  accessible-  conical  intersections  exist.  The  actual 
crossing  seam  occurs  for  geometrical  configurations  in  which  H+  approaches  NO  approximately 
perpendicular  to  the  NO  intemuclear  axis.  These  configurations  are  quite  distinct  from  the 
previously  determined  allowed  crossing  seams  which  correspond  to  collinear  approach.  Thus 
these  results  are  expected  to  have  important  implications  for  the  energy  redistribution  which 
accompanies  charge  transfer.  The  incorporation  of  this  actual  crossing  seam  into  any  set  of 
potential  energy  surfaces  used  to  describe  the  charge  transfer  process  appears  essential.  It  will  be 
quite  interesting  to  see  whether  this  feature  is  to  be  found  in  other  charge  transfer  reactions. 

IV.  FUTURE  DIRECTIONS 
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The  research  completed  during  the  course  of  this  grant  has  suggested  the  need  for 
additional  methods  development  and  applications.  These  future  research  directions  are  indicated 
below.  The  question  of  spin-forbidden  radiationless  decay  in  energetic  species  appears  to  be  a 
fairly  general  one.  The  methodology  developed  as  pan  of  this  grant  for  the  direct  determination  of 
the  minimum  energy  crossing  point  can  provide  considerable  insight  into  the  feasibility  of  a  spin- 
forbidden  process,  particularly  when  this  point  turns  out  to  be  the  barrier  for  this  process. 

However  to  enable  studies  of  the  dynamics  of  these  processes  or  when  the  barrier  does  not 
correspond  to  the  minimum  energy  crossing  point  it  is  necessary  to  determine  more  of  the  crossing 
hypersurface.  To  this  end  a  new  phase  of  program  development  has  been  initiated  with  the  goal  of 
extending  our  capabilities  for  the  direct  determination  of  actual/ailowed/avoided  crossing  hyperlines 
and  h3q)erplanes.  The  algorithms  under  development  will  permit  systematic  determination  of  (i) 
actual/allowed  crossing  seams  passing  through  a  minimum  energy  crossing  point  and  (ii)  avoided 
crossing  seams  parametrized  by  an  arbitrary  set  of  internal  coordinates.  This  program 
development,  to  be  completed  during  our  present  AFOSR  sponsored  research  program,  will 
provide  powerful  new  tools  for  the  study  of  spin-forbidden  and  spin-allowed  electronically 
nonadiabatic  processes  in  polyatomic  systems. 

Recent  theoretical  work  in  other  laboratories  serves  to  motivate  applications  of  the  methods 
developed  in  our  research  program.  Sophisticated  theoretical  methods  have  suggested  the 
existence  of  additional  metastable  geometrical  isomers  of  NO  dimer'^  and  a  metastable  dication 
of  04^^  (isoelectronic  with  N2O2)  has  also  been  predicted.  It  will  be  important  to  learn  whether 
spin-forbidden  radiationless  decay  plays  a  significant  role  in  limiting  the  lifetime  of  these  species. 
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